Mechanism of resonant x-ray magnetic scattering in NiO 



Jun-ichi Igarashi and Manabu Takahashi 
Faculty of Engineering. Gunma University, Kiryu, Gunma 376-8515, Japan 

( ) 

We study the resonant x-ray magnetic scattering (RXMS) around the K edge of Ni in the anti- 
ferromagnet NiO, by treating the 4p states of Ni as a band and the 3d states as locahzed states. We 
propose a mechanism that the 4p states are coupled to the magnetic order through the intra-atomic 
Coulomb interaction between the 4p and the Sd states and through the p-d mixing to the 3d states 
of neighboring Ni atoms. These couplings induce the orbital moment in the 4p band, and thereby 
give rise to the RXMS intensity at the K edge in the dipolar process. It is found that the spin-orbit 
interaction in the 4p band has negligibly small contribution to the RXMS intensity. The present 
model reproduces well the experimental spectra. We also discuss the azimuthal angle dependence 
of the intensity. 
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I. INTRODUCTION 



Resonant x-ray scattering has recently attracted much interest as a useful tool to investigate magnetic and orbiijaiiMj-. 
ders. The resonant enhancement for magnetic Bragg reflections has been observed in transition-metal compoundstlHa^El 
using synchrotron radiation with photon energy around the K edge. Also the resonan|t| enhancement for the prohibited 
Bragg reflection corresponding to the orbital order has been observed in LaMnOs.Q The resonant scattering is de- 
scribed by a second order process that a photon is virtually absorbed by exciting a core electron to unoccupied states 
and then emitted by recombining the excited electron with the core hole. Therefore, for the K edge, the 4p states 
are involved in the dipolar process, while the 3d states are involved in the quadrupolar process. Since the 4p states 
are highly extended in space, they are rather sensitive to electronic structure at neighboring sites. This character is 
now well recognized in the resonant scattering for the orbital order in LaMnOs; at the beginning, proposed was a 
mechanism that the intensity comes from the 4p states modified by the intra-atomic Coulomb interaction|j»ath_the 3c? 
electrons constituting the orbital order but the subsequent studies of treating the Ap states as a bandtO'ElLlllj have 
revealed that the Ap states of Mn is mainly modified by the Jahn- Teller distortion via the oxygen potential on the 
neighboring sites, not by the intra-atomic Coulomb interaction. 

In this paper, we study the mechanism of the resonant x-ray scattering for the magnetic order (RXMS) in NiO. 
We analyze the recent experiment of Neubeck et al.p paying special attentions to the band effects of the Ap states. 
In NiO, Ni atoms form an fee lattice. It is a type II antiferromagnet (AF) with the Neel temperature Tn = 523 K. 
The order parameter is characterized by a wavevector Q, which is directed to one of four body-diagonals in the fee 
lattice. Therefore there exist four domains [K domains). The spin vector lies in the plane perpendicular to the AF 
modulation direction so that there are three possible directions {S domains). For example, the spins are directed to 
(-1,-1, 2), (2, -1, -1), or (-1, 2,-1) in the K domain of Q = (i, \, \). In total, we have 12 domains. In the RXMS 
experiment for NiO, the K domain is specified, but the S domains are mixed. □ 

For generating the RXMS intensity, orbitals have to coupled to the magnetic order; in the dipolar process, the 
Ap states have to couple to the 3c? states. Since such couplings are usually much smaller than those to the Jahn- 
Teller distortion in LaMnOs, the mechanism is much subtler than that for the orbital order in LaMnOa. The spin 
polarization alone in the Ap band cannot give rise to the RXMS intensity; if the amplitudes of exciting electron are 
summed for up and down spins, the total becomes the same at each site in the absence of the orbital moment. This 
leads to the cancellation of the amplitudes between different magnetic sublattices. The spin-orbit interaction in the 
Ap band may break this cancellation,Ba but its [effect will be shown to be negligibly small. 

As discussed in our previous paper for CoOjHl we propose that two couplings are effective to make the Ap states 
couple to the magnetic order: the intra-atomic Coulomb interaction between the Ap and 3c? states, and the mixing 
of the Ap states with the 3c? states of neighboring Ni atoms. To substantiate this observation, we use a tight-binding 
model for the Ap states, treating the 3c? states within a cluster model. _We will show that these couplings actually give 
rise to the RXMS spectra in agreement with the recent experiment A caution is to be made on the intra-atomic 
Coulomb effect; without treating the Ap states as a band, the Coulomb interaction gives rise to too large RXMS 
intensities. Thus the Ap band effect has a role of reducing the effect of the Coulomb interaction. We have also a 
clear peak in the pre- if -edge region, which will be shown to come from the quadrupolar process. Since the 3c? states 
strongly couple to the magnetic order, iaving this peak is not surprising. The present cluster model analysis, a slight 
extension of an atomic model analysis ,□ reproduces well the experimental spectra. 

This paper is organized as follows. In Sec. II, the present model is described, and the electronic structure is calculated 
on the basis of this model. In Sec. Ill, the formulation of the RXMS spectra is developed, and the calculated results 
are discussed in comparison with experiments. Section V is devoted to concluding remarks. In the Appendix, the 
geometrical factors are given in explicit forms. 



For describing the 3c? states mixing to the 2p states of neighboring O sites, we start by the following model 
Hamiltonian: 



II. 3D STATES 






(2.1) 
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Operators dima and pjia- represent the annihilation of an electron with spin a in the 3d state m of the Ni site i 
and of an electron with spin a in the 2p state £ of the O site j, respectively. The first term describes the energy 
of the 3d level. We take account of a point charge crystal-field splitting (lODq) such that Ed{eg) = Ed + QDq, 
Ed{t2g) — Ed — ADq. The second term, H^^_^^, represents the intra-atomic Coulomb interaction between the 3d 
states on Ni sites. Although the explicit form is omitted here, the matrix elements are written in terms of the Slater 
integrals F'^{3d, 3d), F^{3d, 3d), and F'^{3d, 3d). The third term represents the spin-orbit interaction of the 3d states. 
The fourth term represents the energy of the 2p level of O. The relative energy is given by the charge-transfer energy 
A defined by A = Ed - E2p + nU for the d" configuration with U = F°(3d, 3d) - ^F^(3d, 3d) - ^F^{3d, 3d). The 
last term represents the mixing of the Ni 3d states with the O 2p states. The mixing parameters t^^^^'l are defined 

by the Slater-Koster integrals, i£idcr)2p,3(i, (pd7r)2p,3ti.0 Most parameter values arCj-assumed to be consistent with a 
previous cluster-model analysiscJ and our previous electron correlation caiculation.Ej They are listed in Table ^. We 
evaluate the atomic values within the Hartree-Fock (HF) approximation £3 For F°(3d, 3d), we use a rather smaller 
value than the atomic one by considering screening effects in solids, while for F2(3d,3d) and F'^{3d,3d) we use the 
atomic values with multiplying a factor 0.8. 

The ground state is calculated within this Hamiltonian as follows. We assume that each cluster consisting of a Ni 
atom and surrounding six O atoms is subjected to a molecular field coming from neighboring clusters; the exchange 
energy in the mean-field approximation is given by H^x — '^(Q)(S) • S, where S represent the spin operator for the 
3d states, and J(Q) is the Fourier transform of the exchange constant with the AF modulation vector Q. Since the 
Neel temperature Tn is given by Tn = | J(Q)|S'(S' -I- l)/(3fcB) in the mean field approximation, the exchange energy 
is expressed as 

- - -kBT^S,, (2.2) 



2S'(S'+1) 

with the r] axis pointing to one of the directions (—1,-1,2), (2,-1,1), or (-1,2,-1). We put 5=1, |(S)| = 1 for 
Ni atoms, and use the experimental value for Tn. Then we make numerical calculations, assuming that the ground 
state is given by a linear combination of the 3d^ configurations (45 states) and 3d^L configuration (360 states) in each 
cluster, where L indicates a state that a hole is in the 2p states of neighboring O sites. Diagonalizing the Hamiltonian 
matrix (405 x 405 dimensions), we obtain the ground state that the weight of being in the 3d^ configuration is 0.865 
with the orbital moment (L,^) = 0.316. This_value is consistent with the value-iL,^) = 0.32 in the recent polarization 
analysis of nonresonant magnetic-scatteringtj and in the LDA+U calculation.!!!! Note that the calculated value of the 
orbital moment is sensitive to the assumed value of lODq. 

In the quadrupolar process of the RXMS, an electron is excited from the Is state to the 3d states in the intermediate 
states. Such excited states are numerically calculated in a similar way to that for the ground state. We diagonalize 
the Hamiltonian matrix within the space of the ls^3d^ configuration (10 states) and the ls^3d^'^L configuration (36 
states), assuming that an attractive interaction V^d ^ 4 eV works between the 3d electron and the core hole only in 
the core hole site. The energy of the core-hole is adjusted such that a peak in the pre-K-edge region is located at 
huj = 8330 eV in the RXMS spectra (see Sec. IV). 

III. 4P STATES 

The RXMS involves a virtual process of a transition from the Is state to the 4p states at Ni sites in the dipolar 
process, so that the Ap states are to be studied in detail. This section is devoted to the study of the 4p states in 
connection with the K-edge absorption spectra. 

We use a tight-binding model for the 4p states. Although many other states such as 4s state at Ni sites and 3s 
state at O sites are mixed in such a rather-high-energy region, this model makes sense as an effective model for simply 
taking account of the 4p band effects; the 4p states keep their atomic character inside the atomic sphere even in solids. 
As relevant couplings to the magnetic order, we introduce the intra-atomic Coulomb interaction between the 4p and 
the 3d states, and the mixing of the 4p states with the 3d states in neighboring Ni atoms. Thereby a part of the 
Hamiltonian concerning the 4p states is given by 

ina <i,i' > crnn' 

+ K-Zd + E [t'tnf^^lnA'rna + H .c) , (3. 1) 
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where Omcr represents the annihilation of an electron with spin a in the 4p state n on Ni site i. The first and the 
second terms describe the_4p bands. The transfer integral t'in'vn' written in terms of the Slater-Koster parameters, 
{pp(j)ipAp and (pp7r)4p 4p.ll3 We have the 4p band width ~ 20 eV, assuming {ppTr)4p^4p ~ 3.5 eV for the nearest- 
neighbor Ni-Ni distance, and the half of this value for the next-nearest-neighbor distance, with an empirical relation 
(pp7r)4p_4p ~ —(j)pa)4p^4p/A. The value of E4p is adjusted such that the peak in the K-edge absorption is located at 
huj — 8437 eV. The third term represents the spin-orbit interaction in the 4p states. The HF atomic value is used for 
<^4p. The fourth term H^p_^^ represents the intra-atomic Coulomb interaction between the 4p states and the 3c? states. 
The matrix elements are written in terms of the Slater integrals, F°(4p, 3d), F^(4p, 3(i), G^(4p, 3d), and G'^(4p, 3c?). 
By the same reason as in H^^_^^, the HF atomic values are used for F^(4p, 3d), G^(4p, 3d), and G'^(4p, 3d) with 
multiplying a factor 0.8. The value of i^°(4p, 3d) is rather ambiguous, since it can be absorbed into the value of i?4p. 
The last term describes the mixing of the 4p states with the 3d states of neighboring sites. The mixing parameter 
^'in Vm is written in terms of the Slater-Koster parameters, {pda)4p^'id and {pd'n:)4p^'id- We assume (pdi7)4p_3d ~ — 1.0 
eV together with an empirical relation {pdTT)4p^zd ~ — (pdcr)4p^3c(/2. Although this mixing is much smaller than the 
transfer between the Ap states, it is not neglected in the RXMS spectra, since it makes the 4p states couple to the 
magnetic order. The parameter values are listed in Table |. 

We numerically calculate the Ap band as follows. The Hi together with Hq is represented in each Ni site within the 
space of states that one electron is in the Ap states with the 3d electrons in the ground state (6 states) and that a 3d 
electron is added to become the 3d^-configuration (10 states). Such states in different Ni sites are coupled through 



the second and the last terms in Eq. (3.1), so that they form a band. Considering two Ni atoms in the unit cell due to 
the AF order, the Hamiltonian matrix has 32 x 32 dimensions for each wave vector, which is numerically diagonalized. 

In the final state of the K-eAge absorption and the intermediate states of RXMS, there is a hole in the Is state, 
which attracts the Ap electron. Assuming that the core-hole potential works only within the core-hole site, we solve 
the potential problem with respect to a local Green's function defined by 

no„.;0«^(w) - / (r(5j,(<)ao„^(t)aj„,^,6o..(0)))e*'^*di, (3.2) 

where bo^ is the annihilation operator of the Is state with spin a at site 0, (• • •) the average over the ground state, 
and T the time ordering. Taking account of the multiple scattering due to the core-hole potential, we have 

no„.;o„.(a;) = J2 <l;0„"."(^)[^ - ^4pn(°)(a;)]o-\,,,.o„,^,, (3.3) 

with 

(0) , , {Ona\R){R\On'a') 

Here the state |_R) represents the Ap band with energy discussed above, ancLecorc is the core hole energy. The 
life-time broadening width F and the core-hole potential V4p are set to be 1.7 eVjfl and 3 eV, respectively. 

The spectral function of the local Green's function, I(lo) = 'Ylin(y^~^ I ''^^^^^'^'niyfiriai'^) ■, is proportional to the K- 
edge absorp|t|ion spectra and the fluorescence spectra. Figure]^ shows /(w) in comparison with the K edge fluorescence 
experiment D The solid and the broken lines represent the spectra for V^p — and 3 eV, respectively. The attractive 
interaction enhances the intensity in the low-energy side. It consists of two peaks. The peak in low-energy region 
corresponds to the if-edge peak in the experimental spectra. We are mainly interested in the energy region of this 
peak in the RXMS spectra. Th peak in high energy region may correspond to the peak around hio = 8365 eV in the 
experimental spectra, but the calculated peak is too large. This discrepancy may come from using the Ap tight-binding 
model; in reality, it is likely to become much broader through the mixing to other states such as the 4s state of Ni 
and 3s states of O. 



IV. RESONANT X-RAY MAGENTIC SCATTERING 



A. Procedure of calculation 



The magnetic scattering amplitude has been reviewed by several authors We summarize here the expression 

of the amplitude in a form suitable to the present calculation. The scattering geometry is shown in Fig. photon 
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with frequency uj, momentum and polarization /x (= cr or tt) is scattered into the state with momentum ky and 
polarization /z' (= a' or tt'). For the system with two magnetic sublattices A and B in the AF order, the cross section 
for elastic scattering at magnetic Bragg peaks is given by 



da 



(x|T;L^,(G,c.)-r^%,(G,^)| , (4.1) 



with 



+ j;:^^,(G,..) + i;;_^,(G,c.). (4.2) 

where G = kj — kj is t he s cattering vector. 

The first term in e£L ( [4.2[ ) represents a non-resonant term, where m and c are the electron mass and the velocity of 
photon, respectively.o A" and B are given by 

A" = A' - (A' • G)G, A' = -4sin2 e{e' x e), (4.3) 
B = e' X e + (k/ X e')(k/ • e) - (k, x e)(k, • e') - (k/ x e') x (k, x e), (4.4) 

where e and e are the initial and scattered polarizations, and kj, kj, and G are normalized vectors of k^, k/, and 
G. L''(G) and S''(G) stand for the Fourier transform of the orbitai and spin angular momentum densities, which 
values are known from the recent non-resonant RXS experimentMN |L(G)|(|S(G)|) = 0.26(0.82), 0.24(0.60), and 
0.12 (0.20) for G = (i, i, i), (|, |, |), and (|, |, |), respectively. We insert these experimental values into Eq. (^). 
The local moments of orbital and spin are give by L''(G) and S''(G) with G — > 0. The above experimental values 
are consistent with the present ground-state calculation in Sec. II. 

The second term in Eq. ( |4.2| ) represents the resonant contribution from the dipolar process. It is expressed by a 
second-order process that a photon is virtually absorbed by exciting the Is electron to the Ap states and then emitted 
by recombining the excited electron with the Is core hole. The corresponding amplitude is expressed as 



j;_^,(G,c.) = ^(P"^')„AC,(c.)P„^„ (4.5) 



aa' 



with the scattering tensor 



The geometrical factors P^, P'^ are given in explicit forms in the Appendix. In Eq. ( |4.6| ), c^Ag = {Eh ~ Eg)/fi, and 
Xa's stand for the dipole operators, xi — x, X2 = y, and xs — z in the coordinate frame fixed to the crystal axes. 
The ground state \g) has an energy Eg. The intermediate state |A) consists of an excited electron on the 4p states 
and a hole on the Is state with energy E\. In the presence of the core-hole potential, it is convenient to rewrite the 
scattering tensor in terms of the local Green's function defined in Sec. III. We have 

M2^,{iu) ^ AlmulJ2^0aaMaA^), (4.7) 

a 

where site is to belong to sublattice 77, and TilOc « 8350 eV. The is the dipole matrix element defined by 
= (4p|r|ls) = P4p(r)rPis (r)r^dr with Ris{r) and Rip{r) being the radial wavefunctions for the Is and 4p 
states. We evaluate it as Ap — 1.52 x 10~^^ cm within the HE approximationllj in the id^ configuration of a free 
Ni ion. The amplitudes at A and B sublattices are nearly cancelled out, leading to an antisymmetric form for the 
difference: 



AP^(u;) -M^(w) = \ a 5 | , (4.8) 
with a, h, and c being some complex numbers. 
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The third term in Eq. (4.2) represents the resonant contribution from the quadrupolar process. The corresponding 
amphtude is expressed as 



7'' 



(4.9) 



with the scattering tensor, 



^7VH = f^E 



mul,g{g\z-,\K'){K'\z^,\g) 
huj- (Ea' ~ Eg) + iT/2' 



(4.10) 



The geometrical factors, Q^, Q'^ , are given in expHcit forms in the Appendix. In Eq. ( [4.10| ), k is the wavevector 



of the incident (and scattered) photon, ^ 4.2 x 10 cm 



{E^i — Eg)/h, and z^^'s stand for the quadrupole 



operators, zi = (\/3/2)(x^ — y^), Z2 = (l/2)(3z^ — r^), Z3 ~ ^/Syz, Z4 — ^/3zx, and Z5 — VSxy in the coordinate 
frame fixed to the crystal axes. The intermediate states |A') consist of an excited electron on the 3d states and a hole 
on the Is states with energy Ea>, which are evaluated by following the procedure described in Sec II. The scattering 
amplitude NJI^,{uj) contains the square of the quadrupole matrix element Ad = {3d\r^\ls) = R^d{r)r^ Ris{r)r^ Ar , 
with i?3d(r) being the radial wavefunction of a free Ni ion. We evaluate it as Ad = 2.10 x 10^^° cm^ within the HF 
approximation in the id^ configuration. Similar to the dipolar process, the scattering tensors N^[ijj) and N^[ijj) are 
nearly cancelled out, leading to an antisymmetric form for the difference: 



/ d e f g\ 

—d h p q 
— e —h r s 
-f -p -r t 
V -9 -q -s -t Oj 



(4.11) 



with d t being some complex numbers. Note that this difference wouli-jvanish if the spin-orbit interaction in the 
3d states were absent. This is consistent with the prediction by Loversey.Eil 



B. Calculated results 

In this subsection, we discuss the RXMS spectra in comparison with the recent experiment of Neubeck et aZ.I 
In their experiment, the K domain of Q = (^,^,^) is selected, but the contributions from different S domains 
are not separated. Fortunately the relative volumes for three S domains have been determined such that 0.55 {Si), 
0.4(^2), 0.05(^3) for Q = (i,i,i), 0.63(5i), 0.37(^2), 0(^3) for Q = (|,|,i), and 0.26(5i), 0.33(^2), 0.41(^3) for 
Q = (|, |, |). In the following, we average the spectra over three S domains multiplying the weights of these relative 
volumes. 

Figure ^ shows the spectra thus calculated as a function of photon energy for G = (5, 5, 5). The azimuthal angle ip 
is set to be zero such that the scattering plane contains the (1, —1, 0) crystal axis, corresponding to the experimental 
situation. The left panels show the spectra for the a ^ a' channel. Only one peak is found at Uu) — 8330 eV in the 
pre-if-edge region, in agreement with the experiment. Since the dipolar process is forbidden in this polarization, this 
peak arises from the quadrupolar process. A strong anti-resonant dip is seen in the low energy side of the peak, which 
comes from an interference of the non-resonant process with the quadrupolar process. The experimental intensity in 
the high energy side of the peak is much smaller than the calculated intensity. This discrepancy may be related to 
the lack of absorption correction on the experimental spectra, which correction should make the intensity in the high 
energy side the same as that in the low energy side. This observation conversely suggests the size of the absorption 
correction. The right panels show the spectra for the cr — > tt' channel. Similar to the a ^ a' channel, we find a sharp 
peak at fiLo — 8330 eV, which arises from the quadrupolar process. The intensity for the energy far below fiLO = 8330 
eV is much smaller than the value for the a a' channel. In addition, we have several structures in the high-energy 
side of the peak. As plotted by the long-dashed line, the non-resonant process and the quadrupolar process cannot 
give rise to sufficient intensities, indicating that the intensity mainly comes from the dipolar process. We have checked 
the effect of the spin-orbit interaction in the 4p band by calculating the spectra with neglecting that interaction; no 
visible change is found. We have a peak around huj = 8346 eV, in agreement with the experiment. To this peak, 
the intra-atomic Coulomb interaction between the 4p and the 3d states contributes substantially (compare the solid 
line with the dashed one). The experimental intensity, which is smaller than the calculated one, should be enhanced 
by absorption correction. We should emphasize here that the dipolar contribution becomes more than 10 times as 
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large as the present values, topJarge in comparison with the experimental one, if the transfer of the 4p electron to 
neighboring sites is neglected.u Therefore, the Ap band effect is very important for reproducing the experimental 
spectra. A large peak at huj 8365 eV is regarded as an artifact of using the tight-binding model, as mentioned in 
Sec. III. 

Figure || shows the spectra as a function of photon energy, at -0 = for G = (|, |, |). The left panels show the 
spectra for the a —> a' channel. The spectral shape is quite similar to that for G = (i, i, i). The right panels show 
the spectra for the cr — s- tt' channel. The non-resonant intensity is larger than that for G = (i, ^, ^). We find a sharp 
peak at hu = 8330 eV, which arises from the quadrupolar process. In the high-energy side of this peak, we also have 
several structures arising from the dipolar process; a peak around huj = 8346 eV agrees well with the experiment. 

Figure ^ shows the spectra as a function of photon energy for G = (|, |, |). The azimuthal angle is now set 
to be 270°, corresponding to the experimental situation. The right panels show the spectra for the a ^ n' channel. 
For the peak at hui ~ 8330 eV, an anti-resonant dip now moved to the high-energy side. We also find a broad peak 
around huj ~ 8340 eV, owing to the dipolar process. 

Figure ^ shows the spectra as a function of the azimuthal angle ■0 for G = (^i 5)1 comparison with the 
experiment. The lowest panel shows the intensity of the peak at huj = 8346 eV in the main-if-edge region for the 
a ^ tt' channel. We multiply a factor such that the calculated values at V' = becomes close to the experimental 
one. The variation resembles to a dependence of sin^ consistent with the dominant contribution of the dipolar 
process. The top and middle panels show the intensities of the peaks at huj = 8330 eV in the pre-if-edge region for 
a ^ a' and the a ^ -k' channels, respectively. We multiply a same factor in both cases such that the values at ?/; = 
become close to the experimental ones. The variation in the top panel looks simply a 90°-shift of the curve in the 
lowest panel. This is a little surprising, since we expect a more complex shape from the quadrupolar process. The 
curve in the middle panel has an extra peak. This may be related to the fact that the dipolar contribution remains 
in this peak. Such an extra peak is not clear in the experiment. 

V. CONCLUDING REMARKS 

We have studied the RXMS around the K edge in the antifcrromagnet NiO, treating the 4p states as a band and 
the ?>d states as localized states. We have proposed a mechanism that the 4p states are coupled to the magnetic order 
through the intra-atomic Coulomb interaction between the Ap and 3ci states, and through the p-d mixing to the Zd 
states of neighboring Ni atoms. Thereby the orbital moment is induced in the Ap states, giving rise to the RXMS 
intensdty in the dipolar process. The calculated spectra in the present model have reproduced well the experimental 
ones.Q Different from previous analysesJj'Q the effect of the spin-orbit interaction in the Ap band is found negligibly 
small within the present model. This reminds, us of the magnetic circular dichroism (MCD) for the AT-edge absorption 
in the ferromagnetic metals Fe, Co, and Ni.E3 It is known that the MCD intensity is mainly generated by the orbital 
moment induced in the Ap band through the mixing to the Zd states of neighboring atoms, not through the spin-orbit 
interaction in the Ap band.Ell i-i 

We have also obtained a peak in the pre-AT-edge region, in agreement with a previous analysis in the atomic model.u 
This arises from the quadrupole process, showing a contrast to the pre-AT-edge peak in the-eesonant scattering for 
the orbital order in LaMnOs. The latter is predicted to be generated by the dipolar processjlj where the states of p 
symmetry (seen from the Mn atoms) have weights in the pre-AT-edge region by mixing to the 3(i states of neighboring 
Mn atoms. Usually the intensity of the quadrupolar process is much smaller than that of the dipolar process, but 
the dipolar contribution cannot become large in the present case because of a weak coupling of the Ap states to the 
magnetic order. 

Furthermore, we have discussed the azimuthal-angle dependence of the spectra. The curve for the pre-AT-edge peak, 
which comes from the quadrupolar process, looks simply a 90°-shift of the curve for the main-A'-edge peak, which 
comes mainly from the dipolar process. 

Finally, we would like to stress again that the tight-binding model is a first approach for taking account of the Ap 
band effects. More realistic band calculations are desirable for obtaining more quantitative results. 
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APPENDIX A: GEOMETRICAL FACTOR 



We summarize here the expressions of geometrical factors for the dipolar and quadrupolar processes. The detailed 
derivation is given in Ref. [11]. For the scattering vectors G=(i,i,i), (|,|,|), and (|, |, |), the geometrical factor 
P'' {fi = a or tt) for incident photon in the dipolar process is given by 



iP^)3 



(cos f3 cos + sin ijj) / V2, 
(cos f3 cos ip — sin -tjj) / \/2, 
— sin /3 cos V', 

[sin 0(cos /3 sin — cos -tp) + cos 9 sin /3] / V2, 
[sin 0(cos (3 sin 'ip + cos Tp) + cos sin /3] / V2, 
(— sin sin /3 sin ^/^ + cos 6* cos 



(Al) 



(A2) 



where with /3 = arccos y/l/S, and 9 represents the Bragg angle. The azimuthal angle ip ^ corresponds to the 
scattering plane containing the (1,-1,0) crystal axi s. T he ge ome trical factor P'*' (/i' = a' or tt') for the scattered 
photon is given by replacing sin 6* by — sin 6* in Eqs. ( |Al[ ) and (A2). 

Also, the geometrical factor (fj, — a or tt) for incident photon in the quadrupolar process is given by 



(Q")3 



— cos 9 cos /? cos 2^ — sin 9 sin /3 sin -0, 
(\/3/2) cos 6* sin^ /3 sin + {V3/2) sin 9 sin 2/3 cos V', 

— cos 9 sin /3(cos (3 sin 2?/; + cos 2ip) / \f2 

- sin ^(cos 2/3 cos "0 — cos /3 sin 'p)/^/2, 
{Q'')a = -cos6'sin/3(cos/3sin20 - cos2i/')/V2 

- sin 6'(cos 2P cos 0; + cos f3 sin ■ip)/V2, 
cos 6i[l - (1/2) sin^ (3] sin 2tP - (1/2) sin 6* sin 2/3 cos?/', 

— (1/2) sm29 cos /3 sin 20 — cos 26* sin /3 cos?/', 

— ( V3/4) sin 26* sin^ /3 cos 20; - (\/3/2) cos 26* sin 2/3 sin ^0 

- (%/3/2)sin26'[l- (3/2)sin2/3], 

(g'^)3 = (l/2V2)sin26'sin/3(cos/3cos20- sin20) 
+ cos 20(cos 2/3 sin ')p + cos /3 cos -0) / %/2 

- (3/4V2)sin26'sin2/3, 

(Q^)^ = (l/2\/2)sin2e'sin/3(cos/3cos20 + sin20;) 
+ cos 26'(cos 2/3 sin 0) — cos /3 cos -0) / V2 

- (3/4\/2)sin26'sin2/3, 

(g'^)j. = -(l/2)sin26'[l - (l/2)sin2/3]cos20> 

+ (1/2) cos 26* sin 2/3 sin 0-- (3/4) sin 26* sin^ /3. 



(A3) 



(A4) 



The geometrical factor Q''^ {/ i' = cr' or tt') for scattered photon is given by replacing sin 9 and sin 29 by — sin 9 and 
— sin 26*, respectively, in Eqs. (A3) and (A4). 
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TABLE I. Parameter values for the tight-binding model of NiO in units of eV. 



Parameter" Parameter 



6.5 F". 2.5 



dd "-J J^pd 



9.8 F^a 1-9 



72 

dd 

fL 6.1 0.68 

A 6.0 Gpd 0.59 

{pda)2p,3d -1.4 (ppcr)4p,4p 3.5 

(pdn)2p,3d 0.63 {ppn)4p,4p -0.88 

WDq 0.50 (pdcr)4p,3d -1-0 

{pdn)4p,3d 0.50 

Csd 0.096(0.12=) C4p 0.098 



" Values in the ls^3d^ configuration. 

Values in the ls^3d^4p^ configuration. 
" Value in the Is^SdP configuration. 
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FIG. 1. Finite section of a NiO crystal. Only Ni atoms are shown. Wave vector Q = (|, i, i) characterizes an antiferromag- 
netic modulation direction. Three arrows in a plane perpendicular to Q indicate possible spin directions {S domains). Spins 
for solid, grey, and open circles align alternately. 

FIG. 2. Spectral function I{lo) for the local Green's function for the 4p states (lower panel), in comparison with the K-edge 
fluorescence spectra (upper panel).'' Solid and broken lines represent the spectra calculated with the core- hole potential Vip = 
and 3 eV, respectively. 

FIG. 3. Scattering geometry of x ray scattering. Incident photon with wave vector k, and polarization cr or tt is scattered 
into the state with wave vector k/ and polarization a' or tt' at Bragg angle 9. The sample crystal is rotated by azimuthal angle 
V' around the scattering vector G = k/ — kj. 

FIG. 4. Magnetic scattering spectra for G = (i, i, i) at ip = 0°, as a function of photon energy, in comparison with the 
experiment. The upper panels show the experimental spectra without absorption correction.'* The left and right panels show 
the intensity for the a —>■ a' channel, and that for the a —f n' channel, respectively. In the right and lower panel, the dashed 
line represents the intensity by neglecting the intra-atomic Coulomb interaction between the 4p and the 3d states, and the 
long-dashed line represents the contribution from the non-resonant and the quadrupolar processes. 

FIG. 5. Magnetic scattering spectra for G = (f , |, f) at tp — 0° , as a function of photon energy, in comparison with the 
experiment. The upper panels show the experimental spectra without absorption correction.* The left and right panels show 
the intensity for the a a' channel and that for the a —^ iv' channel, respectively. 

FIG. 6. Magnetic scattering spectra for G = (|, |, |) at tp — 270°, as a function of photon energy, in comparison with the 
experiment. The upper panels show the experimental spectra without absorption correction.'' The left and right panels show 
the intensity for the a a' channel and that for the a —>■ n' channel, respectively. 

FIG. 7. Magnetic scattering intensity as a function of azimuthal angle -0, for G = (i, i, i), in comparison with the experi- 
ment. The open circles represent the experimental intensities.'' The top and middle panels are for the peak at hui — 8330 eV 
in the pre-A'-edge region for the a a' and a —^ tt' channels, respectively. The lowest panel is for the peak at Tiw = 8346 eV 
in the main-JsT-edge region for the cr — » tt' channel. 
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